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Optical measurements of the real and imaginary frequency dependent conductivity of uncompen- 
sated n-type silicon are reported. The experiments are done in the quantum limit, hu) > fcsT, 
across a broad doping range on the insulating side of the Metal-Insulator transition (MIT). The 
observed low energy linear frequency dependence shows characteristics consistent with theories of 
a Coulomb glass, but discrepancies exist in the relative magnitudes of the real and imaginary com- 
ponents. At higher energies we observe a crossover to a quadratic frequency dependence that is 
sharper than expected over the entire dopant range. The concentration dependence gives evidence 
that the Coulomb interaction energy is the relevant energy scale that determines this crossover. 



PACS numbers: 72.20.Ee, 71.30.-Fh, 71.45.Gm 

Doped semiconductors close to the MIT are systems 
ideally suited for the study of strong electron-electron 
interactions. Below a critical doping the material is an 
insulator, i.e. at zero temperature there is no DC conduc- 
tivity. The charge carriers are localized in the Anderson 
sense, and screening is reduced from the metallic regime. 

Anderson coined the term Fermi glass to describe a 
non-interacting disordered insulating system whose uni- 
versal properties, independent of system specific details, 
are determined by Fermi statistics alone [1]. Mott's origi- 
nal treatment of such a system [2] , which has a finite den- 
sity of states at the Fermi level, but where localization 
makes an insulator, did not address the issue of electron- 
electron interactions as pointed out by PoUak [3], and 
then Efros and Sklovskii (ES) [4]. A Fermi glass that in- 
cludes interactions between localized electrons has been 
termed a Coulomb glass. Among other things, such a sys- 
tem is typified by a depletion in the single particle density 
of states around the Fermi level which was termed by ES 
the Coulomb gap, A [5]. 

By taking into account the mean Coulomb interac- 
tion between two sites forming a resonant pair U{r^) = 
e^/ei^cj, where — ^[Zn(2/o/?iw)] is the most probable 
hop distance between pairs, /q is the hopping attempt 
rate and £i is the dielectric constant, ES derived the real 
part of the ac conductivity of a Coulomb glass to be: 



CTi = l3e^g^^£,^uj[ln{2Io/riuj)Y[nuj + Uir^)] 



(1) 



Here /3 is a constant of order one, is the non-interacting 
density of states and f is the localization length. This 
formula takes on a different frequency dependence in two 
limits. When the photon energy, Tiuj > U{r^), one re- 
covers the same quadratic frequency dependence that 
Mott derived for a non-interacting Fermi glass. Here the 
Coulomb glass is indistinguishable from the Fermi glass 
in so far as the high frequency limit of the conductivity 
is concerned. In the opposite limit, Hlu < U{r^), the con- 
ductivity of a Coulomb glass will show an approximately 



linear dependence on frequency, plus logarithmic correc- 
tions. The imaginary component of the complex conduc- 
tivity, (72, as predicted by Efros [6], should be identical 
to CTi up to a logarithmic factor. We should note that 
Eq. (1) was derived for the case where huj > A, the 
Coulomb gap width. However a linear dependence (albeit 
with additional logarithmic corrections) and an eventual 
crossover to Mott's non-interacting quadratic law is still 
expected even for the case where huj < A. 

Despite its fundamental importance, there have been 
very few AC conductivity studies that have been done in 
the quantum limit hto > ksT, yet at low enough frequen- 
cies to probe the relevant energy ranges of Eq. (1). Only 
very recently, have measurements been attempted that 
address these issues. M. Lee et al. found that for concen- 
trations close to the MIT the expected linear to quadratic 
crossover occurs, but is much sharper than predicted [7]. 
Other recent work was done on amorphous NbSi, where 
the general systematics of a linear in frequency Coulomb 
glass response was found, but a crossover to the non- 
interacting quadratic dependence was not reached [8] . 

In this Letter we report the first measurements of 
the frequency dependent real and imaginary conductiv- 
ity on the insulating side of the MIT for a doped, crys- 
talline semiconductor. We find a concentration and fre- 
quency dependence of the components consistent with 
predictions of a Coulomb glass, but a discrepancy arises 
in the predicted ratio of their relative magnitudes. At 
higher frequencies a crossover to quadratic Mott-like be- 
havior is observed and due to having measured across 
a broad dopant range, the concentration dependence of 
the crossover energy is shown to be consistent with the 
Coulomb interaction energy U and not the Coulomb gap 
width A. 

Nominally uncompensated n-type silicon samples were 
obtained from Recticon Enterprises Inc. A boule of sili- 
con was grown using the Czochralski method and sliced 
into 1 mm thick slabs. Room temperature resistiv- 
ity was measured using an ADE 6035 resistivity gauge. 
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and the concentration determined using the Thurber 
scale [9]. A number of samples were etched with a 
A%HF + m%HN03 solution; this resulted in no differ- 
ence in the results of frequency dependent conductivity 
measurements. The Si:P samples discussed in this Letter 
span a range from 39% to 69%, stated as a percentage 
ratio of the sample's dopant concentration to the critical 
concentration at the MIT. 
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FIG. 1. T = 2.8 K real part of the frequency dependent 
conductivity versus frequency. The Si:P samples shown are 
at 50% and 69% dopant concentration relative the critical 
concentration, Xc ~ 3.5 x lO^^cm"^. The solid lines are linear 
and quadratic fits to the lower and upper portions of the data 
respectively. The dotted line is a fit to the form of Eq. (1). 

At 35 and 60 GHz the real part of the conductivity 
was evaluated from the measured loss of highly sensi- 
tive resonant cavities via the perturbation method. The 
technique and analysis is well established [10]. In the 
millimeter spectral range, 80 GHz to 1 THz, backward 
wave oscillators were employed as coherent sources in a 
transmission configuration [11]. Fabry-Perot like reso- 
nances in the transmission were analyzed uniquely de- 
termining both components of the complex conductivity. 
Cavity measurements were performed down to 1.2 K and 
in the millimeter spectral range down to 2.8 K. With 
these base temperatures the quantum limit hoj > ksT 
{IK K. 20GHz) of the system was being investigated. 

The real part of the frequency dependent conductiv- 
ity, (Ti, at T = 2.8 K for two samples is shown in Fig. 
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FIG. 2. The real part and the magnitude of the imaginary 
part of the optical conductivity for three different Si:P sam- 
ples on the insulating side of the MIT is shown. The stated 
percentages are the sample's dopant concentration relative 
the critical concentration. The dotted lines are linear best 
fits showing the trend of increasing conductivity as the MIT 
is approached. The error bars shown for tJi are representative 
of those for all the data points. Some of the lower frequency 
points were determined using resonant cavities at a base tem- 
perature of 1.2 K. 

1. This data, representative of all the samples, shows an 
approximately linear dependence at low frequencies and 
then a sharp crossover to an approximately quadratic 
behavior at higher frequencies. This is the qualitatively 
expected behavior from Eq. (1). However, as seen by 
the overlaycd fits, Eq. (1) provides only a rough guide. 
The solid lines are linear and quadratic fits to the low 
frequency and high frequency data respectively. Individ- 
ually these functions fit the data well, however the full 
crossover function does not. The dotted line is a fit using 
the same method as Ref. [7], namely forcing the linear 
portion to pass through the origin as well as the low fre- 
quency data and leaving the pre-factor of the quadratic 
term as a free variable. The crossover between linear 
and quadratic portions is much more abrupt than the ES 
function predicts and is observed over our entire doping 
range as was observed previously in an analogous system, 
Si:B, for samples closer to the MIT [7]. 

Having shown that the frequency dependence is qual- 
itatively consistent with a Coulomb glass, we focus on 
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the low frequency regime. Both the real part and the 
magnitude of the imaginary part, \o2\- of the complex 
conductivity are plotted as a function of frequency in 
Fig. 2. The dotted lines are linear fits and show the rel- 
ative increase in the magnitude of both components as 
the MIT is approached. 

A complex conductivity obeying a power law can be 
expressed in a simple Kramers-Kronig compatible form, 
a\{uj) + i (72(0;) = A{iu})°', in order to have an indepen- 
dent means for determining the exponent. To determine 
the power a one can utilize the ratio of \a2\ versus ai 
(with the frequency as a variable). The power a is given 
by, 



2, -1/^2! 
a = —tan 



(2) 
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FIG. 3. The upper panel shows the ratios of the magnitude 
of the imaginaxy to the real part of the complex conductiv- 
ity for samples of Si:P and amorphous NbSi. The NbSi data 
is adapted from Rcf 8. The middle panel shows the calcu- 
lated powers of a as determined from Eq. 6. The dashed 
line through the NbSi data is a guide to the eye. The bot- 
tom panel shows the divergence of the localization length (us- 
ing lo = lO^^s"'^), and the dotted line is a power law fit of 
(1 — x/xc)~^ , where (3 is found to be 0.95. 

The uppermost panel in Figure 3 shows the ratio of the 

imaginary to the real part of the conductivity for Si:P. 
NbSi is an amorphous insulating glass which is charac- 
terized by a vastly higher density of states (DOS) at Ep 



[8]. Disorder should give a somewhat similar functional 
dependence for the AC conductivity, but with very differ- 
ent energy scales. Data from it is included as another ex- 
ample of AC conductivity in a disordered system and for 
comparison purposes, as Si:P has a very small DOS in its 
impurity band. First we note that this ratio for Si:P re- 
mains large and approximately constant across our range 
of dopant concentrations. From theory, one expects \ai\ 
to be approximately equal to |cr2| to within a factor of 
2-5 (with a reasonable estimate for /q) as predicted by 
Efros [6]. Applied to Si:P, Eq. (1) in the huj < U{r^) 
limit correctly predicts a linear correspondence between 
CTi and I (72 1, but the theory incorrectly predicts the mea- 
sured proportionality by at least a factor of thirty. The 
proportionality is closer for NbSi but has a dependence 
on the doping concentration. This may be related to 
entering the quantum critical (QC) regime as discussed 
below. Here we have used the susceptibility due to the 
interacting electrons themselves Airx = ei — esi, namely 
the full measured dielectric constant, £1, minus the static 
background dielectric constant of the host silicon, esi, to 
determine the magnitude of the imaginary component of 
the conductivity shown in Figure 3. The magnitude of 
the ratio becomes even larger and the discrepancy greater 
if the full dielectric constant Si is considered instead. 

The middle panel in Figure 3 shows the power a as 
determined by Eq. (2). The values for Si:P are approx- 
imately equal to, but slightly less than one, consistent 
with Fig. 2. This indicates that the prefactor of the real 
and imaginary components of the complex conductiv- 
ity have the same concentration dependence. When ap- 
proaching the MIT, the frequency dependence is expected 
to cross over to the QC behavior [12,13], i.e. ui oc u)^^'^ 
in NbSi, when the localization length ^ becomes com- 
parable to the characteristic frequency dependent length 
scale, e.g. the dephasing length, [14]. The crossover 
is not a phase transition and need not be sharp, there- 
fore looking at a fixed window of frequencies, a broad, 
smooth crossover from u) — > cj^/^ would show an aver- 
aged power of the frequency dcpcndcnc;c similar to that 
measured for NbSi shown in the middle panel of Figure 
3. The fact that we see an a w 1 across our whole doping 
range in Si:P, but an a that approaches 0.5 in NbSi in- 
dicates that the critical regime in Si:P is much narrower. 
Simple dimensional arguments [15] give a result similar 
to the non-interacting case [16] that the crossover should 
be inversely proportional to the dopant DOS. The vastly 
smaller (a factor of 10'^) dopant density in Si:P relative 
to NbSi is consistent with a narrower QC regime in Si:P 
as compared to NbSi. 

Because our data spans a large range of concentrations, 
the doping dependence of the crossover energy scale from 
linear to quadratic can be analyzed to see whether its de- 
pendence is consistent with the functional form of other 
energy scales, e.g. the Coulomb interaction energy U or 
the Coulomb gap width A. Fig. 4 shows data of our ex- 



3 



1000 




0.2 

1-x/x^ 



FIG. 4. This work's experimentally determined crossover 
energy scale is represented by the full circles with a fit to a 
function of the relative concentration (1 — x/xa) commonly 
used to normalize data from different materials. 

perimentally determined crossover energy scale plotted 
versus the commonly used fimctional form, 1 — x/xc- M. 
Lee et al. performed AC conductivity measurements on 
a pair of Si:B samples [7], and in said measurements, a 
sharp crossover analogous to our own observations was 
seen. Our data extends the observed range of the sharp 
crossover deep into the insulating regime as M. Lee et al.'s 
original measurements were closer to the critical concen- 
tration. 

Recall that the Coulomb interaction energy between 

two sites forming a resonant pair is U{r^) = e^/eir^, 
which is dependent on concentration via the dielectric 
constant and the localization length dependent most 
probable hop distance. The Coulomb gap width as de- 
termined by ES [4] is A = e^gl^^/el^^. Although M. Lee 
et al. postulated that a smaller Coulomb gap governed 
the crossover, it is not unreasonable to expect that this 
gap scales with the ES single particle one as both are pre- 
sumably caused by the long range Coulomb interaction. 
The magnitude of the gap width is dependent on con- 
centration through the dielectric constant and through 
the DOS term. Thus it is clear that a separate measure- 
ment of the full dielectric constant, £i as a function of 
concentration, as we have done, need be performed in 
order to analyze the crossover energy scale. Although 
there is no theoretical expectation for ujc or A to scale 
as a power law over the whole doping range, we may still 
parameterize these quantities over our doping range as 
power laws. We find that the crossover frequency lUc is 
proportional to (l — x/xc)^-^^, and the full dielectric con- 
stant ei oc (1 — x/xc)~^''^- Using these values, it seems 
highly unlikely that the hojc = A, as the required DOS's 
concentration dependence for this to be valid seems im- 
probable. 

We do find strong support that the crossover energy is 



determined by the Coulomb interaction energy though. 
By setting the measured crossover energy scale equal to 
the Coulomb interaction energy we are able to determine 
both the magnitude of the localization length and the 
exponent with which it diverges as a function of concen- 
tration. Using an appropriate pre-factor for the over- 
lap integral [17], Iq = lO^'^.s^^ in the most probable hop 
distance term, r^, we find a localization length depen- 
dence as plotted in the bottom panel of Fig. 3. We find 
that the localization length exponent is close to unity, 
the value originally predicted by McMillan in his scaling 
theory of the MIT [18], and the magnitude of the local- 
ization length is reasonable. Both these results strongly 
point towards the Coulomb interaction energy as being 
the energy scale at which the observed frequency depen- 
dent crossover from ES to Mott-like hopping conduction 
occurs. 

In summary, we have observed behavior consistent 

with a Coulomb glass across our entire range of dop- 
ing concentrations in Si:P for the frequency and con- 
centration dependence of the real and imaginary compo- 
nents of the conductivity. It is consistent with theoretical 
predictions except in predicting the relative magnitudes 
of (Ti and C72 and the sharper than expected crossover 
from linear to quadratic frequency dependence. From 
the expected form of U{x) and the reasonable behav- 
ior of the extracted localization length, it seems likely 
that the crossover is governed by the Coulomb interac- 
tion strength of a resonant pair and not the Coulomb 
gap. 
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